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dr Abstract 

'• 

^—j ' The dependence of cosmological and inflationary parameters on time during 

the last 60 e-folds in inflation is investigated using a slow-roll inflation model. 

1 ^ i ' The time dependence of the inflaton field is calculated for the case of chaotic 

inflation by a numerical method rather than the familiar approximations, ft is 
04 ■ found that the Hubble parameter and the spectral index decrease in the last 60 

e-folds. The dependence of the power spectrum of the curvature perturbation 
on the size perturbation k is calculated numerically, and it is shown that the 
, overall constant value of the power spectrum thus determined differs from that 

t-H ' given by the familiar method. However, the fc-dependent spectra of the nu- 

, merical calculation, the Bessel function approximation, and the familiar Taylor 

— ' expansion are all found to be consistent. 

oo : 

O ■ PACS number:98.80Cq 

> : 

X ' 1 Introduction 
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Inflation is an important concept in cosmology, and is supported by recent 
satellite-based measurements. The present authors have studied the effect of 
the length of inflation and pre-inflation on the power spectrum and the angular 
power spectrum [1], and it has been shown that the suppression of the spectrum 
at Z = 2 as indicated by Wilkinson Microwave Anisotropy Probe (WMAP) data 
[2] may be explained to a certain extent by the finite length of inflation for 
inflation of 50-60 e-folds [3]. The length of inflation is generally believed to 
be very long, but the physically interesting period of inflation is that in the 
last 60 e-folds. Although a correct inflationary potential based on superstring 
or supergravity theory has yet to be established, it could be assumed that the 
potential term in this latter period of inflation can be expressed in a simple form, 
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such as a small-field or large-field model, or a hybrid model. In the present 
paper, the time-dependent behavior of the inflaton field during the last 100 
e-folds in inflation is investigated using a simple inflation model. Specifically, 
numerical calculations are performed using a chaotic inflation model (4> 2 model) , 
and the effects of the time dependence of the inflaton field on cosmological and 
inflationary parameters such as slow-roll parameters, the Hubble parameter, and 
the spectral index are investigated. The /c-dependence of the power spectrum 
of the curvature perturbation has been investigated using Taylor expansion or 
the numerical method [4-7]. Here, using the derived parameters, i.e., the slow- 
rolled parameters, Hubble parameter and the inflaton field which have the time 
dependence, the power spectrum is calculated with the numerical method and 
the Bessel approximation, and it is compared with the power spectra derived 
from the familiar methods. 



2 Cosmological and inflationary parameters 

The dependence of the cosmological and inflationary parameters on time is 
investigated using a slow-roll inflation model. Assuming a spatially flat universe 
and the energy density of the universe to be dominated by the inflaton field, the 
Einstein field equations can be written as 

and 

4> + 3H4> + V'(4>) = 0, (2) 

where <j> is the inflaton field, V((f>) is the inflaton potential, m is the Planck 
mass, and H — a/a with a being a scale factor. The variable N defined by 
dN — Hdt is introduced to represent time with respect to the number of e-folds 
of inflation. Equation (1) can then be written using the relation <j> = Hd<p/dN 
as 

2 = ^V{<t>) 

3m 2 (l -4ir/3m 2 (d(f>/dN) 2 )' y> 

The Hubble parameter is usually considered to be a constant, i.e., H 2 = 8tt / '3m 2 V '(</)) 
. In the generalization, however, with <j> being dependent on time, the Hubble 
parameter also includes a time dependence (from cq. (3)). Using the variable 
N, eq. (2) can be written as 

d 2 cf> 1 dH d<t> d<j) V'U) 

— — H — + 3—— H — =0 (4) 

dN 2 H dN dN dN H 2 { ' 

Equation (4) can be further written in terms of <f> and the derivate of <j> with 
respect to N. To solve eq. (4) numerically, the slow-roll inflation model is 
adopted, which is a chaotic inflation model with potential term given by V{4>) = 
M 4 /2((f>/m) 2 . Using an initial condition of N = no, 



<t>{no) = J^ m ' ( 5 ) 

V 27T 
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d<p(n ) = m 

dN ~ y/8n(-no) ' 

The numerical solutions for eq. (4) for n values of —100 to —60 do not reveal 
any appreciable dependence on n . The behavior of <f> in the case of no = —100 
is shown in Figure 1. The inflaton field can be seen to slowly decrease. The 
behavior of the Hubble parameter according to eq. (3), which includes <j>, is 
shown in Figure 2. The Hubble parameter becomes considerably smaller as N 
approaches zero. The value H 2 (N) / H 2 (-60) is 1.33 for N = -80, 0.835 for 
N = -50, 0.669 for N = -40, and 0.504 for N = -30. For slow-roll inflation, 
the following parameters are employed [4]: 

( ~ 6 2 ( 2 +V) ~ ^ { H{4>) ] ' (7) 
s _m 2 H»(c» 

t _ m 4 H'WH»'(<t>) 



s 16tt 2 {H(<j>)) 2 ' 

Other slow-roll parameters (ey, r]v, £v) can be written in terms of the slow-roll 
parameters e, S, and ^ to the first order in slow roll: e = ev, 3 = Vv — £v, 
and C = - 3ey?7y + 3e^, where e v = m 2 /IQtt{V /V) 2 , J]y = m 2 /8tt(V" /V), 
and C = m 4 /64ir 2 {V'V"' /V 2 ). Figure 3 shows the behavior of e. The value 
of e can be seen to remain small (e(— 60) = 0.0083, e(— 50) = 0.0099, and 
e(— 40) = 0.0124) for N smaller than —10. The value of the important parameter 
5 is zero in the first-order approximation using this cf) 2 model, and the values of 
5 and £ are very small. The time dependence of the spectral index n s can also 
be calculated in this manner. The spectral index is given by 

n s = 1 + 26 - 4e (10) 

and the time dependence of n s is shown in Figure 4. The value of n s decreases 
as N approaches zero: n s (-80) = 0.981, n s (-60) = 0.975, n s (-50) = 0.970, 
and n s (-40) = 0.963. 



3 Scalar perturbations 

The dependence of the curvature perturbation on the parameter k is investigated 
assuming a spatially flat Friedman-Robertson- Walker (FRW) universe with a 
background spectrum given by metric perturbations. The line element for the 
background and perturbations is generally expressed as [8] 

ds 2 = a 2 {T){{l+2A)dT 2 -2d l Bdx l dT-[{l-2^)5 l3 +2d i d J E+h t3 ]dx l dx ] } 1 (11) 

where r is the conformal time, the functions A, B, *f>, and E represent the scalar 
perturbations, and hij represents tensor perturbations. The density perturba- 
tions in terms of the intrinsic curvature perturbation of comoving hypersurfaces 
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is given by 1Z = — — (H/<j))5<f>, where <j> is the inflaton field, d<p is the fluc- 
tuation of the inflaton field, and 1Z is the curvature perturbation. Overdots 
represent derivatives with respect to time t, and primes represent derivatives 
with respect to the conformal time r. Introducing the gauge-invariant potential 
u = o,(t)(54>+ {4>/H)^) allows the action for scalar perturbations to be written 
as [9] 

S=\J dTd*x{(^-f - (Vu) 2 + ^u 2 }, (12) 

where Z = atp/H and u = —Z1Z. The field u{r,x) is expressed using annihila- 
tion and creation operators as 

u(t, x) = ^-L^ J d 3 k{u k (r)a k + 4(r)ai k } e - jkx , (13) 

and the field equation for Uk{r) is derived as 

d 2 Uk „, 1 d 2 Z , „ .„ 

^ + ^-z^ ttfc = ' (14) 

where the solution to Uk satisfies the normalization condition Ukdu* k / 'dr—u^duk/ 'dr = 
i. Using the slow-roll parameters, (d 2 Z/dr 2 )/Z is written exactly as 

where 

v = (l + e-'±S + e 2 -2e5+^ + ^). (16) 
A partial integration allows r to be expanded as follows [10]. 

1 2e(e - S) f e(-2S 2 + 6(9 - 4e)e - 6e 2 + 3e 3 - g(l - e)) dN 

aH(l-e) aH(l-e) + J (1 - e) 2 off' 

(17) 

Here, using the approximate relation a 2 H 2 = (1 + 2e(e — S)) 2 /t 2 (1 — e) 2 , eq. 
(14) can be rewritten as 

If it is assumed that the slow-roll parameters can be constant, the solution for 
eq. (18) can be written using the Hankel function Hj^\—kr) as 

U k (r) = J^ e *(^V2) 1 r/2 ( _ fcr) l/2 ff (l) ( _ fcr)j (19) 



where \i = \j2v(\ + 2e(e - <5)) 2 /(l - e) 2 + 1/4. This solution is fixed such that 
as fcr — > — oo, Mfe(r) approach plane waves. The numerical solution to eq. (14) 
can then be considered. Equation (14) can be written in terms of N as 

d 2 u k 1 daHduk . k 2 

dN 2 - + ^H^NdN + ( JaW = °- (20) 



4 



In order to calculate eq. (20) numerically, it is assumed that the present-day 
size perturbation k = 0.002(1/Mpc) exceeds the Hubble radius in inflation at 
the time of N = -60. The scale factor is thus written as a = 0.002/iJ e Ar+6 °, 
where Ho is the value of H at N = — 60. As the Hubble parameter and v have 
both been calculated as functions of N, eq. (20) can be calculated numerically 
with an initial value of Uk derived from eq. (19). Note that the desired solution 
of Uk is that at the point aH = k. In Figure 5, the absolute values of Uk derived 
by numerical calculation of eq. (20) at aH = k is shown as a function of k. The 
absolute value of Uk (eq.(19)) obtained using the Bessel approximation, adopting 
the value of \i as the numerically calculated value, is also shown. The value of 
fj, changes slightly according to N: /i(— 80) — 1.513, /z(— 60) = 1.517, /i(— 50) = 
1.520, At(-40) = 1.525, and /x(-30) = 1.534. The results shown in Figure 5 do 
not suggest any appreciable difference between the numerical calculation and 
the Bessel approximations, indicating that the Bessel approximation, in which 
\x exhibits k dependence, can be considered satisfactory. 

The power spectrum of the scalar perturbations P-r is defined as follows [4] . 

2tt 2 

<7e k (7 7 ),7er(7 7 )>=^ r p^ 3 (k-i). (21) 

Here, TZk{v) is the Fourier series of the curvature perturbation 1Z. The power 
spectrum P-r is then written as [4] 

fl.-(£)l*l", (22) 

where Z — a<p/H. The power spectrum is estimated by four different methods; 
the numerical method, the Bessel approximation using eq. (19), and Taylor 
expansion with and without running as two familiar methods. The numerical 
calculation of the power spectrum is performed by the present method for k = 
aH. In the case of the Bessel approximation using the values of /z, d<j)/dN, and 
H derived from the numerical calculation (having k dependence), the power 
spectrum (22) can thus be written as 

_ e (-fcr)(gW(-fcr)) 2 
n ~8n a 2 (#/rfiV)2 |fe=aH ' [Z6) 

where the value of Pr is obtained at k = aH (i.e. fcr = — (l + 2e(e— — e))). 
There is very little difference between the results of the numerical calculation 
and the Bessel approximation (see Table 1), again demonstrating that the Bessel 
approximation is satisfactory. In the first of the familiar cases, the power spec- 
trum is written as 

"-^'•jw^-'ffl^swi-" < 24 > 

In the alternate case, the power spectrum is given by [5] 

P n = (1 - (2C + l)e + C8f-^— 2 \ k=aH , (25) 
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where C = — 2 + ln2 + 7 = —0.73, with 7 being the Euler constant. The value 
of the power spectrum (P-ji) given by eq. (24) is 3049.07(M 4 /to 4 ), and that by 
eq. (25) is 3100.34(M 4 /m 4 ) assuming the value of \i at N = -60. Using the 
same parameter values, the numerical calculation by eq. (22) and the Bessel 
approximation (eq. (23)) both afford a value of 6252.43(M 4 /m 4 ). There is thus 
a difference of approximately a factor of two between the present numerical and 
Bessel calculations and the familiar cases (eqs.(24) and (25)). This discrepancy 
can be attributed to the difference in the Hankel function used for deriving 
the expressions. The familiar expressions are derived from the asymptotic form 
(k/aH — ► 0) of the Hankel function, and are estimated at k/aH = 1[4], whereas 
the present Bessel calculation uses the correct form. The difference between 
the asymptotic and correct forms of the Hankel function equates to a factor of 
approximately \[2 at k/aH — 1, giving rise to a change in the normalization of 
the power spectrum. 

As the value of the spectral index differs from unity, as indicated by ex- 
perimental results such as the WMAP data [2, 11], the power spectrum can be 
inferred to have k dependence. The power spectrum is usually expanded at a 
pivot scale ho, i.e. [6], 

Pn(k) = P K (fco)(^)" s(fco) - 1+(1/2)Qln(fe/fco) , (26) 
fc 

where a — dn s /dlnk\ k=ko . Two cases are considered here; a — (no running of 
the spectral index), and a 7^ (running). Table 1 lists the results for the four 
calculations at various values of k, and Figure 6 compares the numerically cal- 
culated power spectrum with that determined by Taylor expansion with power 
spectrum running, where the over-all factor Pjz(ko) is fixed as 6252.43(M 4 /m 4 ). 
These results confirm that the power spectrum given by Taylor expansion with 
running exhibits very similar behavior to that obtained by the present numerical 
calculation, and that there is very little difference among the power spectrum 
values for the numerical, Bessel, and Taylor expansion with running cases. The 
largest difference, 0.05% occurs at k = 1 (1/Mpc), while the Taylor expansion 
without spectral running differs by up to 1%. 

4 Discussion and summary 

The time dependence of the cosmological and inflationary parameters and the k 
dependence of the power spectrum of the curvature perturbation during the last 
100 e-folds in inflation were investigated numerically using a slow-roll inflation 
model. Using a time dependence for the inflaton field calculated for the case of 
chaotic inflation ((j) 2 model), it was shown that both the Hubble parameter and 
spectral index decrease with time (see Figures 2 and 4). The time dependence 
of the slow-roll parameters was also derived. It was shown that an approxi- 
mation of the Hankel function obtained by changing the slow-roll parameters 
is closely consistent with the numerical calculations, and that there exists an 
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approximately two-fold difference in normalization between the present meth- 
ods and the usual methods due to the asymptotic form of the Hankcl function 
in the usual treatments. In the fc-dependent spectra, the numerical and Bessel 
calculations are very similar to the Taylor expansion with spectral running, al- 
though the Taylor expansion without running differs by up to 1% at k = 1 
(1/Mpc). These results indicate that the running of the spectral index should 
be considered in the derivation of an inflation model. Although calculations 
were performed in the present study for chaotic inflation with potential of (j> 2 , 
this approach can be similarly applied to any arbitrary potential of inflation. It 
is inferred that similar behavior as derived in the present case will be obtained 
for other slow-roll models. Detailed analyses using other slow-roll models will 
be performed as part of future work. 
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Table 1: Power spectrum values (P-r.[M 4 /to 4 ]) and ratio to numerical calcula- 
tion 



fc(l/Mpc) 


Numerical 


Bcssel 


Taylor(running) 


Taylor (no running) 


0.002 


6252.43 


6252.43 


6252.43 


6252.43 


0.01 


5923.03 


5923.39 


5921.54 


5925.80 


(ratio) 


(1) 


(1.00006) 


(0.999749) 


(1.00047) 


0.1 


5467.71 


5468.01 


5464.68 


5487.96 


(ratio) 


(1) 


(1.00006) 


(0.999446) 


(1.0037) 


1.0 


5030.45 


5030.75 


5028.23 


5028.23 


(ratio) 


(1) 


(1.00006) 


(0.999559) 


(1.01034) 



0/m 









3.5 
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\ 1.5 



Figure 1: Variation in 4>/m (m is Planck mass) with respect to N. 
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Figure 2: Variation in the Hubble parameter (H 2 (N)/H 2 (— 60)) with respect 
to N. 



e 




I N 

-80 -60 -40 -20 

Figure 3: Variation in e with respect to N. 
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n s 









0.975 








95 








\ 0.925 


-80 


-60 


-40 


~ \ 0.875 
\ 0.85 
01.825 



N 



Figure 4: Variation in n s with respect to N. 



|u k | 



0.2 0.4 0.6 0.8 



1 



k (1/Mpc) 



Figure 5: Variation in Uk with respect to k. Dotted line denotes the numerical 
calculation, and solid red line denotes the Bessel approximation using the nu- 
merically calculated value of \i. 
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Figure 6: Variation in Pn with respect to k. Dotted line denotes the numerical 
calculation, and the red solid line denotes the Taylor approximation with spec- 
tral index running for the same value of Pn(ko)- 
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